Cloud information content analysis of multi-angular measurements in the oxygen A-band: application to 3MI and MSPI

The vertical distribution of cloud cover has a significant impact on a large number of meteorological and climatic processes. Cloud top altitude and cloud geometrical thickness are then essential. Previous studies established the possibility of retrieving those parameters from multi-angular oxygen A-band measurements. Here we perform a study and comparison of the performances of future instruments. The 3MI (Multi-angle, Multi-channel and Multi-polarization Imager) instrument developed by EUMETSAT, which is an extension of the POLDER/PARASOL instrument, and MSPI (Multi-angles Spectro-Polarimetric Imager) develoloped by NASA's Jet Propulsion Laboratory will measure total and polarized light reflected by the Earth's atmosphere–surface system in several spectral bands (from UV to SWIR) and several viewing geometries. Those instruments should provide opportunities to observe the links between the cloud structures and the anisotropy of the reflected solar radiation into space. Specific algorithms will need be developed in order to take advantage of the new capabilities of this instrument. However, prior to this effort, we need to understand, through a theoretical Shannon information content analysis, the limits and advantages of these new instruments for retrieving liquid and ice cloud properties, and especially, in this study, the amount of information coming from the A-Band channel on the cloud top altitude (CTOP) and geometrical thickness (CGT). We compare the information content of 3MI A-Band in two configurations and that of MSPI. Quantitative information content estimates show that the retrieval of CTOP with a high accuracy is possible in almost all cases investigated. The retrieval of CGT seems less easy but possible for optically thick clouds above a black surface, at least when CGT > 1–2 km

Benchmarking clear-sky reflectances

Accurate calculations of shortwave reflectances in clear-sky aerosol-laden atmospheres are necessary for various applications in atmospheric sciences. However, computational cost becomes increasingly important for some applications such as data assimilation of top-of-atmosphere reflectances in models of atmospheric composition. This study aims to provide a benchmark that can help in assessing these two requirements in combination. We describe a protocol and input data for 44 080 cases involving various solar and viewing geometries, four different surfaces (one oceanic bidirectional reflectance function and three albedo values for a Lambertian surface), eight aerosol optical depths, five wavelengths, and four aerosol types. We first consider two models relying on the discrete ordinate method: VLIDORT (in vector and scalar configurations) and DISORT (scalar configuration only). We use VLIDORT in its vector configuration as a reference model and quantify the loss of accuracy due to (i) neglecting the effect of polarization in DISORT and VLIDORT (scalar) models and (ii) decreasing the number of streams in DISORT. We further test two other models: the 6SV2 model, relying on the successive orders of scattering method, and Forward-Lobe Two-Stream Radiance Model (FLOTSAM), a new model under development by two of the authors. Typical mean fractional errors of 2.8 % and 2.4 % for 6SV2 and FLOTSAM are found, respectively. Computational cost depends on the input parameters but also on the code implementation and application as some models solve the radiative transfer equations for a range of geometries while others do not. All necessary input and output data are provided as a Supplement as a potential resource for interested developers and users of radiative transfer models

Liquid cloud optical property retrieval and associated uncertainties using multi-angular and bispectral measurements of the airborne radiometer OSIRIS

In remote sensing applications, clouds are generally characterized by two properties: cloud optical thickness (COT) and effective radius of water–ice particles (Reff), as well as additionally by geometric properties when specific information is available. Most of the current operational passive remote sensing algorithms use a mono-angular bispectral method to retrieve COT and Reff. They are based on pre-computed lookup tables while assuming a homogeneous plane-parallel cloud layer. In this work, we use the formalism of the optimal estimation method, applied to airborne near-infrared high-resolution multi-angular measurements, to retrieve COT and Reff as well as the corresponding uncertainties related to the measurement errors, the non-retrieved parameters, and the cloud model assumptions. The measurements used were acquired by the airborne radiometer OSIRIS (Observing System Including PolaRization in the Solar Infrared Spectrum), developed by the Laboratoire d'Optique Atmosphérique. It provides multi-angular measurements at a resolution of tens of meters, which is very suitable for refining our knowledge of cloud properties and their high spatial variability. OSIRIS is based on the POLDER (POlarization and Directionality of the Earth's Reflectances) concept as a prototype of the future 3MI (Multi-viewing Multi-channel Multi-polarization Imager) planned to be launched on the EUMETSAT-ESA MetOp-SG platform in 2024. The approach used allows the exploitation of all the angular information available for each pixel to overcome the radiance angular effects. More consistent cloud properties with lower uncertainty compared to operational mono-directional retrieval methods (traditional bispectral method) are then obtained. The framework of the optimal estimation method also provides the possibility to estimate uncertainties of different sources. Three types of errors were evaluated: (1) errors related to measurement uncertainties, which reach 6 % and 12 % for COT and Reff, respectively, (2) errors related to an incorrect estimation of the ancillary data that remain below 0.5 %, and (3) errors related to the simplified cloud physical model assuming independent pixel approximation. We show that not considering the in-cloud heterogeneous vertical profiles and the 3D radiative transfer effects leads to an average uncertainty of 5 % and 4 % for COT and 13 % and 9 % for Reff.</p

Ice particle habit and surface roughness derived from PARASOL polarization measurements

Ice clouds are an important element in the radiative balance of the earth's climate system, but their microphysical and optical properties still are not well constrained, especially ice particle habit and the degree of particle surface roughness. In situ observations have revealed common ice particle habits and evidence for surface roughness, but these observations are limited. An alternative is to infer the ice particle shape and surface roughness from satellite observations of polarized reflectivity since they are sensitive to both particle shape and degree of surface roughness. In this study an adding–doubling radiative transfer code is used to simulate polarized reflectivity for nine different ice habits and one habit mixture, along with 17 distinct levels of the surface roughness. A lookup table (LUT) is constructed from the simulation results and used to infer shape and surface roughness from PARASOL satellite polarized reflectivity data over the ocean. Globally, the retrievals yield a compact aggregate of columns as the most commonly retrieved ice habit. Analysis of PARASOL data from the tropics results in slightly more aggregates than in midlatitude or polar regions. Some level of surface roughness is inferred in nearly 70% of PARASOL data, with mean and median roughness near σ = 0.2 and 0.15, respectively. Tropical region analyses have 20% more pixels retrieved with particle surface roughness than in midlatitude or polar regions. The global asymmetry parameter inferred at a wavelength of 0.865 μm has a mean value of 0.77 and a median value of 0.75

An algorithm to retrieve ice water content profiles in cirrus clouds from the synergy of ground-based lidar and thermal infrared radiometer measurements

The algorithm presented in this paper was developed to retrieve ice water content (IWC) profiles in cirrus clouds. It is based on optimal estimation theory and combines ground-based visible lidar and thermal infrared (TIR) radiometer measurements in a common retrieval framework in order to retrieve profiles of IWC together with a correction factor for the backscatter intensity of cirrus cloud particles. As a first step, we introduce a method to retrieve extinction and IWC profiles in cirrus clouds from the lidar measurements alone and demonstrate the shortcomings of this approach due to the backscatter-to-extinction ambiguity. As a second step, we show that TIR radiances constrain the backscattering of the ice crystals at the visible lidar wavelength by constraining the ice water path (IWP) and hence the IWC, which is linked to the optical properties of the ice crystals via a realistic bulk ice microphysical model. The scattering phase function obtained from the microphysical model is flat around the backscatter direction (i.e., there is no backscatter peak). We show that using this flat backscattering phase function to define the backscatter-to-extinction ratio of the ice crystals in the retrievals with the lidar-only algorithm results in an overestimation of the IWC, which is inconsistent with the TIR radiometer measurements. Hence, a synergy algorithm was developed that combines the attenuated backscatter profiles measured by the lidar and the measurements of TIR radiances in a common optimal estimation framework to retrieve the IWC profile together with a correction factor for the phase function of the bulk ice crystals in the backscattering direction. We show that this approach yields consistent lidar and TIR results. The resulting lidar ratios for cirrus clouds are found to be consistent with previous independent studies.</p

Influence of ice particle model on satellite ice cloud retrieval: lessons learned from MODIS and POLDER cloud product comparison

The influence is investigated of the assumed ice particle microphysical and optical model on inferring ice cloud optical thickness (τ) from satellite measurements of the Earth's reflected shortwave radiance. Ice cloud τ are inferred, and subsequently compared, using products from MODIS (MODerate resolution Imaging Spectroradiometer) and POLDER (POLarization and Directionality of the Earth's Reflectances). POLDER τ values are found to be substantially smaller than those from collocated MODIS data. It is shown that this difference is caused primarily by the use of different ice particle bulk scattering models in the two retrievals, and more specifically, the scattering phase function. Furthermore, the influence of the ice particle model on the derivation of ice cloud radiative forcing (CRF) from satellite retrievals is studied. Three sets of shortwave CRF are calculated using different combinations of the retrieval and associated ice particle models. It is shown that the uncertainty associated with an ice particle model may lead to two types of errors in estimating CRF from satellite retrievals. One stems from the retrieval itself and the other is due to the optical properties, such as the asymmetry factor, used for CRF calculations. Although a comparison of the CRFs reveals that these two types of errors tend to cancel each other, significant differences are still found between the three CRFs, which indicates that the ice particle model affects not only optical thickness retrievals but also CRF calculations. In addition to CRF, the effect of the ice particle model on the derivation of seasonal variation of τ from satellite measurements is discussed. It is shown that optical thickness retrievals based on the same MODIS observations, but derived using different assumptions of the ice particle model, can be substantially different. These differences can be divided into two parts. The first-order difference is mainly caused by the differences in the asymmetry factor. The second-order difference is related to seasonal changes in the sampled scattering angles and therefore dependent on the sun-satellite viewing geometry. Because of this second-order difference, the use of different ice particle models may lead to a different understanding of the seasonal variation of τ

From Architectured Materials to Large-Scale Additive Manufacturing

The classical material-by-design approach has been extensively perfected by materials scientists, while engineers have been optimising structures geometrically for centuries. The purpose of architectured materials is to build bridges across themicroscale ofmaterials and themacroscale of engineering structures, to put some geometry in the microstructure. This is a paradigm shift. Materials cannot be considered monolithic anymore. Any set of materials functions, even antagonistic ones, can be envisaged in the future. In this paper, we intend to demonstrate the pertinence of computation for developing architectured materials, and the not-so-incidental outcome which led us to developing large-scale additive manufacturing for architectural applications

On the reflection and polarisation properties of ice cloud

International audienceA number of cirrus ice crystal scattering models are tested using measurements of total reflectance and polarised reflectance obtained from the space-based polarisation and directionality of Earth's reflectances (POLDER) instrument. In this paper, 1 day of global POLDER data is utilised taken from the 25 June 2003 to test the assumed ice crystal models. The POLDER instrument is able to test the validity of various ice crystal models since it can measure the total reflectance and polarised reflectance at up to 14 different viewing directions almost simultaneously between the scattering angles of about 60-180°. It is found that ice crystal models that are randomised (in this case the randomisation element is through distortion) from some pristine ice crystal geometry best fit simultaneous measurements of total and polarised reflectance. The optimal distortion parameter that best describes the POLDER measurements is found to be 0.40, which has been applied to a randomly oriented six-branched bullet-rosette and randomly oriented chain-like aggregate. Moreover, distorted ice crystals that have undergone significant distortions beyond 0.40 may fit the total reflectance measurements but not the polarisation measurements. Therefore, total reflectance measurements by themselves do not provide sufficient information to constrain assumed complex/distorted ice crystal models

A self-consistent scattering model for cirrus. I: The solar region

International audienceIn this paper a self-consistent scattering model for cirrus is presented. The model consists of an ensemble of ice crystals where the smallest ice crystal is represented by a single hexagonal ice column. As the overall ice crystal size increases, the ice crystals become progressively more complex by arbitrarily attaching other hexagonal elements until a chain-like ice crystal is formed, this representing the largest ice crystal in the ensemble. The ensemble consists of six ice crystal members whose aspect ratios (ratios of the major-to-minor axes of the circumscribed ellipse) are allowed to vary between unity and 1.84 for the smallest and largest ice crystal, respectively. The ensemble model's prediction of parameters fundamental to solar radiative transfer through cirrus such as ice water content and the volume extinction coefficient is tested using in situ based data obtained from the midlatitudes and Tropics. It is found that the ensemble model is able to generally predict the ice water content and extinction measurements within a factor of two. Moreover, the ensemble model's prediction of cirrus spherical albedo and polarized reflection are tested against a space-based instrument using one day of global measurements. The space-based instrument is able to sample the scattering phase function between the scattering angles of approximately 60° and 180° , and a total of 37 581 satellite pixels were used in the present analysis covering latitude bands between 43.75°S and 76.58°N. It is found that the ensemble model phase function is well able to minimize significantly differences between satellite-based measurements of spherical albedo and the ensemble model's prediction of spherical albedo. The satellite-based measurements of polarized reflection are found to be reasonably described by more simple members of the ensemble. The ensemble model presented in this paper should find wide applicability to the remote sensing of cirrus as well as more fundamental solar radiative transfer calculations through cirrus, and improved solar optical properties for climate and Numerical Weather Prediction models

Three-dimensional polarized Monte Carlo atmospheric radiative transfer model (3DMCPOL): 3D effects on polarized visible reflectances of a cirrus cloud

International audienc